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Recent advances in lattice and continuum QCD have given us new insights into quarko- 
nium physics. These set of lectures are intend for the uninitiated. We first give a 
' physical picture of quarkonium and describe the hybrids states established in lattice 

QCD. Then we give an unorthodox presentation of Non-Relativistic QCD (NRQCD) 
MH ■ including a novel method for the application of spin-symmetries. Finally we describe 

the prototypical application of NRQCD: cancellation of infrared divergences in decays 
of P-wave quarkonia. 
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1 Introduction 



> 

• 1.1 About these lectures 

(N 

I have put together a little material on the theory of quarkonia that I hope will 
be a good introduction to the subject for novices, be them graduate students or 
researchers outside this field. I must warn the readers that I give a personal view of 
the subject. There are many excellent, more orthodox reviews available. However, 
I have combined material in a way that I think brings out the physics more clearly. 
For example, although my aim is to present the advances in NRQCD of the last few 
O |' years, at the expense of the whole subject of quarkonium production I have used 

much space and effort in giving a physical picture of quarkonium and hybrids based 
on the rather good recent calculations of inter-quark potentials in lattice QCD. I also 
have injected here and there little tricks that help me understand or calculate. For 
example, I have given explicitly a foolproof method of calculating spin-symmetry 
relations using the Wigner-Eckart theorem by means of a trace formula, akin to 
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C3 ' that used in HQET. 



There are three parts to these lectures. First we discuss the physical picture. To 
this end we briefly use the bag-model picture to motivate a non-trivial inter-quark 
potential, linear at long distances. Then the role of octets in the description creeps 
in through the discussion of hybrids. Hybrids are the particle physics analogs of 
molecules, so it should be no surprise that a Born-Oppenheimer approximation can 
be used to describe them (the BO is reviewed). This immediately suggests we study 
QCD in the Non-Relativistic limit, or 'NRQCD', which is the subject of the middle 
part of these lectures. The last part describes the application of NRQCD to the 
decays of S and P-wave quarkonia. 

Since these lectures are intended for learning students, I have toiled to include 
many exercises. They are included as separate paragraphs, and most often the 
answer is given. 

I have left many loose ends untied. This is not only because I have been careless. 
There are many questions that are unresolved, and I hope the alert reader will spot 
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them and hopefully solve them. In the conclusions I attempt to list some of these 
questions, but offer little in the way of solutions. 




Fig. 1. Static quarks at small separation. The shaded region represents 
the non-perturbative QCD vacuum. The 'bag' of perturbative vacuum 
has radius given in Eq. (0) 



1.2 Quarkonium: A Physical Picture 

Consider a quark Q and an antiquark Q pinned down a distance r from each other. 
Let us begin our considerations by assuming that this distance r is much smaller 
than the typical size of a hadron Aq^ d ~ 1 Fermi. See Fig. |l| Far away from the 
pair, at a distance R 3> r, the chromo-electric field is that of a dipole, 

r cos 6 

Ee~g{r)— £3-, (1) 

where g(r) is the coupling constant of QCD at distance scales of order r. This de- 
scription in terms of a perturbative vacuum cannot be correct at very large distances 
from the dipole. When the field strength drops to a critical value Eg tCri t ~ Aq CD 
the vacuum quickly turns non-perturbative. Thus the die/ale can be described as a 
perturbative "bag" inside the non-perturbative vacuum.Su'u The radius Rf, of this 
bag is thus determined 
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Fig. 2. At short distances the inter-quark force is Coulomb-like, and 
arises from this one gluon exchange diagram. 



We will be interested in dynamical bound states, with quarks free to move. For 
this we need an understanding of the force between static quarks as a function of 
distance r. For r <C Rt the interaction is well approximated as a Coulomb force. It 
arises from perturbative single gluon exchange; see Fig. |[ The energy of the state 
is the energy of the Coulomb dipole, ~ 1/r. At the opposite extreme, shown in 
Fig. |3[ when r 3> Rb non-perturbative effects become important. The perturbative 
vacuum is squashed out into a cigar shape; this is the "string limit" . The bag carries 
volume energy density £ . If the cross sectional area of the cigar in the string limit 
is A, the energy of the state is £ ■ (Ar), ie, it increases linearly with r. The linear 
slope of this potential energy is the "string tension" £A. a 




Fig. 3. For large inter-quark separation r the perturbative vacuum 'bag' 
takes on a cigar shape, of fixed cross sectional area A. The energy of the 
configuration is proportional to the volume of the bag, Ar, giving rise 
to a linear inter-quark potential. 



tn fact the cross sectional area of the string solution is not an independent parameter, A ~ 
g 2 /£, so the string tension ~ yj g 2 £. 



1.3 Digression: The strength of the Coulomb Gluon 

In the perturbative vacuum the inter-quark force is dominated by single gluon ex- 
change, as in the following figure: 



3,P 



i,p' =p + q 



4 Y T^T t a k 



l,-p 



k, -p> 



Fig. 4. The gluon exchange gives rise to a Coulomb interaction. 



The momentum dependence gives rise to a Coulomb potential (the Fourier transform 
of l/q 2 is proportional to 1/r). Here we will focus on the color factor, ^o^ij-^lfe' 
The matrices T° are the generators of SU(3) in the fundamental representation, 
that is, Gell-Mann matrices. The indices i, j, k, I refer to the external quarks, as 
in Fig. ||. 

It is now easy to see that the strength of the interaction depends on the relative 
color of the QQ state. Since Q is in the 3 (fundamental representation) of SU(3), 
and 

3x3 = 108, 

the QQ force is in either of these two channels. It will be instructive to carry 
the calculation for the more general case of SU(N), so the representation 8 is 
understood as the adjoint representation of SU(N). To hnd the 1 and 8 components 
of the force we prepare the initial QQ state to be purely 1 or 8: 

(QjQl)i = Sji(QQ) (3) 
(QjQih^QjQi-jjSjtiQQ). (4) 

Contracting this with 

a 

we get 

(QjQ^J^TtjTll 

a 

(Q j QihJ2 T tj T it 

a 



N 2 -l 



2N 
1 

'm 



(QiQk)l 



(QiQ k ) 8 - 



(5) 
(6) 
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We see that the force in the color octet channel is of opposite sign (in fact, 
repulsive) and weaker than in the singlet channel by a factor of 1/(N 2 — 1). 

In atomic physics the L ■ S interaction is computed most easily via the well 
known formula 

(2L-S) =j(j + l)-l(l+l)-s(s + l). 
The reader may verify the analogous formula 

£ T a T a ) R = C(R) - 2C(fund) 

a 

where C (R) is the Casimir invariant of the representation R, 

^ T a T a = c ^ lm 

a 

For SU(S) if we normalize C(R) so that C(8) = 3, then C{R) = 4/3, 10/3 and 6, 
for R — 3, 6 and 10, respectively. 




Fig. 5. One gluon exchange gives rise to gluon-quark force. The force is 
attractive in the 3 and 6 channels, but repulsive in the 15 channel. 



One can compute the force between a gluon and a quark similarly; see Fig. 
The available channels are given by 

3x8 = 15060 3. 

The force is attractive in the 3 channel, 

C(3) - (7(8) - C(3) = -3, 

and in the 6 channel, 

C(6~) - C(8) - C(3) = -1, 
but repulsive in the 15 channel. 

1.4 Hybrids: Octet Quarkonia 

As we have seen a quark- antiquark pair in a color octet configuration repel each 
other, and a quark-gluon pair attract each other in either the 3 or 6 configurations. 
A state of quarks in an octet configuration can bind to a gluon, provided the gluon 
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attraction to both quarks overcomes the repulsion between quarks. States like this in 
which the glue state plays a particle-like role in binding are called "hybrids" , because 
they are a hybrid of a glueball and a pure quarkonium state. Lattice calculations El 
convincingly indicate that QCD predicts such states, but experimentally no such 
state has been demonstrated. The state is very non-perturbative in nature, but can 
be empirically understood in terms of a bag picture in the string limit with excited 
chromo-electric and magnetic cavity modes in the bag. 

Hybrids have a useful analog in atomic physics: the Hi molecule. There two 
heavy protons repel each other, but the intervening electronic cloud attracts them 
both and overcomes the repulsion. The physics of the molecule is well described by 
the Born-Oppenheimer (BO) approximation, which we now review. 

1.5 Quick Review: Born-Oppenheimer Approximation 
Consider i the molecular Hamiltonians 

H = T e + T N + V 

where the kinetic energies for electrons and nuclei are, respectively, 

d 2 „ ^ 1 d 2 



T, 



1 * 7 J J 



and the potential energy V = V(Xj,Xi) is a function of all particle locations, in 
general. 

In the BO approximation one considers first the simpler Hamiltonian 

H {0) =T e + V 

obtained formally by taking the infinite nuclear mass limit. Since [Xj, H^] = we 
can simultaneously diagonalize Xj and H (°) . Set the values of Xj to Xj and notice 
that X' appear as parameters of the Hamiltonian for the electrons. We include a 
label X' in eigenstates and eigenvalues, 

H^\nX') = W n {X')\nX'), 

to remind us that they carry implicit dependence on X' . In the \xX) representation, 

(xX\nX') = </) n {x, X')S(X - X') 

this is just the Shrodinger equation for the electrons in the field of fixed point 
charges 

[T e + V(x,X')](f> n (x,X') = W n (X')(f> n (x,X'). 

Next we let the nuclei move, slowly. We replace above X' — > X'(t). If the 
motion is slow enough we can solve the new Schrodinger equation 

[T e + V(x,X'(t))]4> n = W n 4> n 
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in the adiabatic approximation: b 



4> n = cj) n (x,X'(t)) 

W n = W n (X'(t)) . 



(up to a phase) 



The condition for the validity of the adiabatic approximation is that the probability 
P n of electrons jumping out of state \nX') be small: 



Pn 



iAW 



«1, 



(7) 



where AW stands for the electronic energy spacings. 

To estimate P n , let a be the molecular size, that is, the typical separation 
between nuclei. By the uncertainty principle the electronic momentum is p e ~ 1/a 
and thus AW ~ E e ~ p 2 /m ~ 1/ma 2 . For the numerator we argue as follows. To 
remove one atom from the molecule we need to move that nucleus by SX' ~ a. At 
this point the new state is essentially orthogonal to the original one, 



for AX' 



dx(j> n (x,X' + AX')(j> n (x,X') w 
or, since </>„ is normalized, 

r f>^ 
dx a 



dX> 



Now, denote the slow nuclear velocity by Vj — dX'Jdt. We have, roughly, 



dt 



-V„ 



1 t_N_ 

aV M 



where the rms velocity is V TI 
Thus 



E,-(^/) and the kinetic energy tjv = E, Wi^f) 



t N /Ma 2 _ / Tn \ % 
^(1/ma 2 ) ~ VmJ £V 



Now, tjv has contributions from rotational and vibrational modes. In most cases 
the latter is dominant, and we write -E v ib = Since the molecule is bound by the 
electronic cloud, the vibrational potential energy V N = \Mu?X 2 at X - a should 
be comparable to the electronic energy, ^Miv 2 a 2 



tN ~ -Evib 



1 



1 1 

2 ma 2 1 

1 



\jMm a 2 



Using this in P n 



m m 
" ~ M V M' 



(8) 



'The phase, here ignored, can play an important rolel 
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Exercise: Show that rotational modes give P n ~ (m/M) 2 . 

In the adiabatic approximation n is a good quantum number (the coupling 
between electronic levels is neglected). Eigenvectors of H are linear combinations 
of \nX') with n fixed: 

J \nX')ij;{X')dX' , 

or, in the \xX) representation, 

<Z> n (x,X)=<t> n (x,X)iP(X) 

The eigenfunctions of H are found by expanding in terms of these. Let's use these 
as trial wavefunction in the variational method. The energy functional is 

FY* ) = (*nl g l*"> = 
K U> (<&n|#„) 

where 

H n ^{X)= J cp n {x,X)[H^ n {x,X)i;{X)]dx, 
and the condition 8E($ n ) — is equivalent to 

H n ijj = Erp 

This is nothing but the Schroedinger equation for the nuclei in the background 
Coulomb field of the electrons in the n-state. 
Exercise: Show that 

H n = T N + W^ + W n 

where 

so the nucleus moves in a potential W n (X) + W' n {X). Show W' n ~ 1/Ma 2 and is 
therefore a small correction. 

1.6 Born-Oppenheimer in Quarkonium and the Relevant Distance Scales 

It seems obvious that the Born-Oppenheimer approximation should be adequate in 
the description of hybrids. There we would choose the quark and antiquark to be 
the slow degrees of freedom while the glue would be the fast degrees of freedom. One 
would guess that the probability P g of Eq. ((?]), that the glue stays in a fixed state, 
is approximated by Eq. (||), P g ~ (Aqcd/-^q) 3 ^ 2 , so for large Mq the adiabatic 
approximation will be very good. 

It is not so obvious that the color singlet states can also be treated in this 
approximation. This is a subtle issue which I won't discuss, save for the following 
guess based on atomic physics analogs. Since the (QQ)i acts like a bound state 
in a mutually attractive Coulomb field, one could guess that the average kinetic 
energy tgg ~ 1/Ma 2 ~ "Rydberg", where a ~ (a s M) _1 is the "Bohr Radius", 
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and if the gluonic energy (or "mass") is fixed E g ~ m g ~ Aqcd, then P g ~ 
(rn g /M)(tQQ/E g ) ~ (AQCD/Af)( a sAf/AQCD) ~ af- Therefore the approximation 
will work to the extent that a 2 is small. 

The scales that appear in this discussion play an important role throughout, so 
we pause to revisit them. The overall mass of the system is dominated by M; the 
size of the system is its Bohr radius, a, with a -1 ~ Ma a ~ Mv; the "Rydberg" is 
the gross measure of energy spacings and is the (inverse of the) primary time scale 
in the system, ~ Mv 2 ; and the strong interactions become relevant at distances 
greater than 1/Aqcd- Just as in the Hydrogen atom, we we have used v ~ a s as 
the velocity of the quarks. So we have the hierarchy of scales 

M >M(j» Mv 2 > A QCD . 

Whether this is realized in nature for charm or beauty is a subject much debate, 
so we will develop the theory hoping that it will be appropriate for at least one of 
these systems. Notice also that theoretically one could have such large masses that 
even Mv n 3> Aqcd for some n > 3. Such scales dictate the physics of, eg, hyperfine 
splittings. In such a world one would be able to compute reliably even such small 
effects. 

It is now clear how our discussion of the static quarks in a perturbative bag 
plays a role: it just corresponds to the first step in the BO approximation. As ar- 
gued above, the inter-quark potential is Coulombic at short distances and increases 
linearly with separation at long distances. In the BO approximation this potential 
is used in a Schroedinger equation for the quarks to describe the bound state. 

1.7 Towards NRQ CD 

Now that we have decided to apply the Born-Oppenheimer approximation to heavy 
quarkonia we find difficulties implementing this program. Consider the Lagrangian 

£ = -0(j 7 • D - M)ij) 

where the covariant derivative is = — i^A^. Dirac's Hamiltonian density 
follows 

H = -tp(-ij ■ d + -A Q j° - -7 • A - M)ip. 
c c 

We immediately face two difficulties. In the standard application of Born-Oppenheimer 
we need to write a Hamiltonian that includes a kinetic energy term for the slow 
degrees of freedom. Our Hamiltonian has no obvious kinetic energy term, and it 
does not account separately for the Q and Q variables. 

The solution to this problem is to work with a non-relativistic version of H, 
by doing a l/c expansion. This immediately leads to what is known as "Non- 
Relativistic QCD" or NRQCD. 

Now, the l/c expansion in QED is old hat. You can find detailed discus- 
sions in famous texts.Q In fact, the Born-Oppenheimer approximation and the non- 
relativistic expansion-applied to QQ states in QCD is also old hat, first considered 
two decades agoJaHMil It will become clear what new developments have re-ignited 
interest in this field. 
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2 NRQCD 



I will present NRQCD as an expansion in a parameter, namely 1/c. The unattractive 
alternative is to expand in a kinematic variable, v Te \ — the relative velocity of the 
QQ pair in a slowly moving CM frame. 

2. 1 Notation 
Our metric is 

f^" = diag(+ ). 

The Dirac 4-spinor is broken up into two 2-spinors 

-CO- < 9 » 

The Dirac gamma matrices satisfy 

and are given, in the Dirac representation, by 



7° 



1 / a* 

-1 I ' 1 1 -cr 4 



where the 2x2 Pauli matrices a 1 satisfy {cr\ o^} = 26^ and a 3 = diag(l, —1). We 



will also use 

1 

1 



■ 1 2 S 

75 = il 7 7 7 = 



and the invariant four index tensor has 



Finally, 

is, in the Dirac representation, 



,0123 _ , _ i 
e — — £oi23 — 



o* _ -„ ; o * _ f *£T* \ « _ _ ij-fc / a fc 



ff =^7 =( ^ ^=-.e- Q ^ 



2.2 Classical NRQCD 

In the Dirac representation the two 2-spinors in Eq. ([)[) roughly correspond to 
the particle and anti-particle components of the 4-spinor. We'd like to make this 
explicit by rewriting the Dirac lagrangian in terms of variables appropriate for a 
non-relativistic approximation about a frame of reference. The appropriate frame 
is, of course, that for which the particles are at, or almost at, rest. To this end shift 
the field so that the trivial (but dominant) time evolution factor is taken out, 

* = e~ iMc2 **, 
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in the Dirac Equation 

[ij ■ D + Mc(7° -1)]* = 



and write out the components * 



\%D t -iff -3 \ U 

iff-D -\(2Mc 2 + iD t )J \X 



0. 



The field \ is suppressed by a power of l/c so one may solve this equation recursively. 
Formally, we write 

y = = iff ■ Dtp 

2Mc+\iD t 

and use the field 

to describe the quark Q being nearly on-shell, in the rest frame of the QQ system, 
in the l/c expansion. 

Exercise: Carry out a similar expansion, starting from 

and show how to describe Q in terms of <f in the 1 /c expansion. 

We can now easily compute the lagrangian that describes the quark state, 

cC = c*[h • D + Mc(7° - 1)]* 

= ib^iDtib — ib^iff ■ D iff ■ Dib 

r w r 2Mc+HD t 

c 1 

This is not yet written as an expansion in 1 /c. To this end we expand 

(2Mc+-iD t )- 1 = (2Mc)" 1 (l - W t /2Mc 2 + ■ ■ •) 
c 



and use 



t''t< /),/), = (S ij +ie ijk a k )D i D j 



= D 2 + ^ff-B 
2c 

to write a 1/c-expanded lagrangian 

cc = ^(iD t ^Df w + ^_fff ■ m + ■ ■ ■ 

Exercise: The first term in the ellipsis is ff ■ DD t ff ■ D. Consider the piece 
ff ■ Da i [D t ,Di\. Show it gives D ■ E + ia ■ D x E. 

Exercise: Let A = ff-D. Write AD t A = \A[D U A]-\[D U A] + \D t A 2 + \A 2 D t . 
As in the previous exercise, the first two terms give no time derivatives D t acting 
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on ijj. We'd like to have time derivatives appear in the lagrangian only as ip'Dti/j. 
Show that this can be accomplished, to this order, by redefining ip — > (1 — ^A 2 )i(j. 

Exercise: Determine the 1/c-expanded lagrangian for the antiquark Q. 

Putting the result of these exercises together, and incorporating the lagrangian 
for the fast degrees of freedom we have the following lagrangian for NRQCD: E3 

C = \{d i Al--d t A a i - 9 -f abc A\A c f 
2 c c 

— — dj^i — —fabcA^A'j) 2 

2Mc V 



8M 3 c 2 

cpg 
~8M 2 c 2 
cs9 



ip^(D ■ E — E ■ D)ip 
ip^a- (Dx E-E x D)%1) 



8M 2 c 2 

+0{l/c 3 ) (10) 

Here we have rescaled ip by \fc and we have generalized the coefficient of the fourth, 
sixth and seventh lines with hindsight (although, for now, cp = Cd = cs = 1). We 
have omitted the antiquark terms, which are left to the student as an exercise. The 
terms in this lagrangian have well known physical interpretation. The fourth line 
is the chromo-magnetic moment interaction, the fifth line is the first relativistic 
correction, the sixth line is the Darwin term and the seventh is the spin orbit 
coupling. 

It would be desirable to have a lowest order NRQCD lagrangian that would be 
completely independent of 1/c. One would then use the 1/c = spectrum as a solid 
basis for perturbation theory in the small parameter 1/c, much as it is done in the 
1/M expansion for B-mesons. If this cannot be done, that is, if our lowest order 
lagrangian contains explicitly the parameter 1/c, then the lowest order states will 
depend implicitly and non-trivially on 1/c, and an explicit expansion in 1 /c will be 
impossible. 

What are the leading terms in the 1/c expansion? To this end one is tempted 
to omit from the lagrangian above all terms with powers of 1/c. This is almost 
correct except for one subtlety. The second term on the first line of ([h]) must be 
kept. Recall, for a Hamiltonian we need the generalized momentum 

dC 1 

7Ti = ■ r = -^(d t Ai - cdiA Q ) 

d(d t Ai) c 2 

so that the Hamiltonian starts at order c 2 : 

H = ndtAi -£= icV + ■ ■ • 
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Also, in the third line, the covariant derivative contains hidden dependence on 1/c 
which should be dropped. Hence, we have the leading order lagrangian: 

2.3 Interlude: Momentum Space Derivation 

It is instructive to derive the effective lagrangian of NRQCD by analyzing Green 
functions (Feynman diagrams) in momentum space. This may be less elegant that 
the derivation above, but it has several advantages. First, it is more intuitive, since 
most of us have spent a great deal of time computing Feynman diagrams. Second, 
it is a better starting point towards including quantum effects. Third, it's easier. 




p — k p 
Fig. 6. Feynman diagram giving the leading contribution to Compton 
scattering of a gluon off a quark. 



Consider Compton scattering of a gluon off a quark; see Fig. |g. The momentum 
p of the internal propagator will be conveniently taken to have components 

p= (Mc+ -E,p), 
c 

This choice is the right starting point for a 1/c expansion; we have written the 
time-like component as an energy, and shifted it by the dominant term Mc. Now, 
we expand the intermediate propagator about 1/c = 0; we assume the momentum 
p "scales" with c as given, ie, we take the large c limit keeping M, E and p fixed: 

1 7 • p + Mc _ Mc(l + 70) + \E^ - 7 • p 
cp 2 - {Mc) 2 ~ (Mc + \E) 2 -p 2 ~ M 2 c 2 
= Mc(l + 70) + ■ ■ ■ 
2ME - p 2 + ■ ■ ■ 

2 J E-p 2 /2M V 

The factor (1 + 7o)/2 in last line projects out the ip component of the 4-spinor, 
and the second factor, 1/(E — p 2 /2M), is the non-relativistic propagator. This is 
precisely the propagator for ip that follows from C ■ 
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For this expansion to work we need to insist that the momenta of any quark 
scales with c just like p does, so both p and p— k scale at most as p = (Mc + ^E,p), 
and therefore k must scale at most as k — (-ko,k). 




Fig. 7. Quark-gluon vertex. In NRQCD the vertex is spin independent. 



Consider next interactions. Since the propagators on either side of the vertex in 
Fig. 0have a projector operator P + , where P± = (l±7o)/2, = P±, P + +P- = 1, 
the vertex can be simplified, 

P+(ig s T a ^)P + = ig s T a 5^ la P+ = ig s T a P+S^. 

Therefore only Aq couples to the quark, and the coupling is a term in the lagrangian 

5£ = ^g 8 T a i>Al 

The propagator and interaction of the quark are given therefore by the Lagrangian 

c = ^m + ^, 

which is precisely the leading term in the quark lagrangian in the previous section. 

We can now try to construct an effective field theory in this limit. I like ±q 
think of an effective field theory as a factorization theorem for Green functions: 13 

g( Pq = (Mc+\E,p),k a = [\k°,k) \ =A(c)G((E : p),(k : k);A?j+O(±). 

There is much explaining to do. The left hand side is the Green function for the 
full theory (QCD), with the kinematic variables as chosen, ie, in the scaling regime 
discussed earlier. The variables pq and k g are generic for the 4-momenta of quarks 
and gluons, of which there may be several. The expression at right shows how the 
full Green function factorizes into a term, A(c), that contains all the dependence on 
the small parameter, and a term that contains the dynamical information, G. What 
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makes the effective theory useful is that G can be computed as the corresponding 
Green function of a theory based on an "effective" lagrangian. 



k 2 




p + ki p + ki + k 2 



Fig. 8. That a diagram like this, with an arbitrary number of gluons, 
is reproduced correctly by NRQCD follows from repeated application of 
the analysis for the propagator and vertex diagrams. 



In fact, we have just proved this identity at tree level for the class of Green func- 
tions which contain one quark line and any number of external gluons, as displayed 
in Fig. ||. 

Let's explore what interesting things may happen at higher order in the loop 
expansion and in higher order in 1/c, as well as for other classes of Green functions. 



2-4 Loops 

Let us take a convergent example, like the Feynman diagram of Fig. ^. This avoids 
unnecessary confusion from ultraviolet divergences. We have seen that the prop- 
agators of QCD go over into those of NRQCD if the loop momentum scales as 
k = (-k ,k). However this means that the gluon propagator becomes, in leading 
order in 1/c, 



1 1.2 _ U2 U2 

c 2 «-q ft, ft, 

This is problematic, since it truncates the theory in a drastic way: it suppresses all 
retardation effects. However, it seems that if we keep the t-A:q in the propagator 
we fail to organize the expansion in powers of 1/c. 
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Fig. 9. This one loop diagram is convergent. One must show that 
NRQCD reproduces QCD's non-analytic behavior at small momenta. 



The solution to this puzzle is to measure internal gluon momentum in energy 
units. That is, we take the gluon momentum to scale as k — {^ko,^k) where ko 
and k have units of energy. This is not inconsistent with our previous scaling rules, 
it is just more restrictive. The propagator is then 



U2 _ U2 



and the loop integral is over 



dk 



d 3 k 



This modification to our scaling rules introduces a surprise. The internal quark 
propagators are now 



E + k - 



(p+^fc) 2 

2M 



E + fcn - tStt 



.1, 



2M 



In NRQCD this further expansion of the quark propagator is called a "multipolc 
expansion", because it is just what we would do if the gluon were real and we 
expanded in inverse powers of the large wavelength A ~ l/|fc|. We emphasize that 
the previous, tree level, success remains valid. 

It is worth dwelling a bit on why we insist in retaining the ko in the gluon 
propagator, rather than following our nose through the 1/c expansion and using 
the propagator in ([Tl|). The problem is that this amount to making static our fast 
moving modes in the Born-Oppenheimer approximation. There are, for sure, static 
components to the interaction — in Coulomb gauge the potential Aq is instanta- 
neous. But using (fyl) removes all retardation effects. 



2.5 Other Green Functions 

Since we are primarily interested in bound states of heavy quark and antiquarks, 
we should study almost forward Q — Q scattering amplitudes; see Fig. ml. 



1G 



(Mc+±E,p) 



{Mc+\E',p' 



U=(HE'-E),p'-p) 



Fig. 10. One gluon exchange. The gluon is harder than the internal 
(loop) gluons in NRQCD. 



The external initial and final quark momenta are chosen according to our scaling 
rules, p = (Mc + -E,p) and p' = (Mc + -E' \p'). This implies the momentum 
transfer q = p 1 — p = (-(£" — E),p' — p), carried by the gluon has a spatial 
component that scales as c° rather than 1/c. Is this a problem? 




Fig. 11. Local interaction reproducing the hard gluon exchange of Fig. llOl 



Our rules state that the internal integration variables for gluon momenta should 
have units of energy, and therefore come accompanied by appropriate factors of 
1/c. But in the diagram of Fig. [l^ the momenta in the gluon propagator is fixed 
externally and scales as c°, while the energy scales as 1/c. So the amplitude is 




This can be reproduced by our effective theory as a spatially non-local 4-quark 
vertex. Since the Fourier transform of l/\q\ 2 is l/\x\, the effective lagrangian should 
be augmented by terms that are schematically of the form 

tpm^x) I d 3 y-J—xh(t,y)- 

J \x-y\ 

This interaction is responsible for, among others, the Coulomb-like potential be- 
tween Q and Q. 

Note that the interaction is local in time, as it must for an appropriate Hamil- 
tonian formulation of the problem. 
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Exercise: Compute explicitly the Coulomb term that needs to be added to the 
effective theory. 



.6 The NRQCD Lagrangian Revisited 



After our exploration in momentum space we return to a derivation of the effective 
lagrangian directly in configuration space. Our aim is to understand in terms of a 
Lagrangian formulation the origin of the new rules that we derived by considering 
Feynamn diagrams. From here on we will adopt Coulomb gauge V ■ A = 0. As we 
saw earlier the lowest order lagrangian density 

A) = l&AZ - - c d t Alf - i(d^) 2 + tf{iD t + |^ 



still contains c dependence. Defining 







Ai(y= ~,t) = y/cAi(x,t) 
c 

and writing the lagrangian in terms of this new variable, we obtain 



L = / d 3 y 



i(d I,) 2 -i(cVy 2 



^{iD t + —)^ + -{d i A f 



2M 



This is independent of c and can be used as a starting point for an expansion in 
1/c. 

In particular: 

• The states of Lq are c-independcnt 

• These states can be used to formulate perturbation theory in 1 j c 

• All 1/c corrections can be expressed as matrix elements of operators between 
these states. That is, we have explicit power counting. The engineering di- 
mensions of an operator can be determined, and compensated for with the 
appropriate inverse powers of mass, M, and speed of light, c. 

Consider the 1/c corrections. The terms of order 1/c in the Lagrangian density 
of Eq. © are 



L c -i = / d 3 x tfi{x,t) 



9 j. (I* A — 



CF9 

2Mc 5 / 2 



tp(x, t), 



where Bi = eijk-^jA]~(y,t). As we see, in these interaction terms the fields have 
dependence both on x and on y = ^x. Changing variables to y does not eliminate 
the dependence on c. Instead we write all field dependence in terms of x and 
expand in a. Taylor series the terms that depend on ^x. This is just the multipole 
expansion: 

d 



L c -i M p = I d 3 x tp f (x,t) 



\ Mc 3 / 2 
cf9 
+ 2Mc 5 / 2 



Ai(0,t) + -^-(djAi(0,t) + ■■■) 



dx l 



u ■ B(0, t) 



(12) 
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Although this expression seems unfamiliar, it is the configuration space version of 
the multipole expansion performed earlier in momentum space were it appeared in 
a natural and intuitive form. 

2. 7 Technical Note On Explicit Power Counting 

The multipole expansion makes power counting straightforward, and this is why 
we have introduced it here. Let us see how it works in more detail in a somewhat 
non-trivial example. 




Fig. 12. Contribution to the quark self-energy from two insertions of the 
operator if)' a ■ Bip, represented by a shaded circle. Without a multipole 
expansion this graph gives contributions of lower order in 1/c than the 
order of the operators inserted. 



Consider a contribution to the quark self-energy from two insertions of the 
operator ijj^a ■ Bip. This is one of the terms of order 1/c in the lagrangian density 
in Eq. ©. 

If we use this operator in the lagrangian as it stands in Eq. (|Io|), that is, before 
multipole expanding, the one loop contribution to the dimensionally regularized 
two-point function is 

T ( 2 ) _ c 2 F C(R)g 2 r jPk (k x g) ■ (k x g) 

~^^ F ~ J W° P + ie)(E + k -^+ ie) ' 

With D = 4 — 2e, the divergent pieces as e — ► are, 

2 (e - + AMc 2 + -^— . 

V 2M J 3 2Af 

Other schemes yield analogous results. 

As an exercise the reader should compute the loop integral with a cut-off 
|fc| < A. Clearly he or she will find a term of order c F c ^ a ° A 2 . This result, and 
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others like it, has led many (including me) to make the ridiculous statement that 
the effective theory must be defined with a low cut-off, say, A ~ Rydberg ~ 1/Ma 2 . 

This is, of course, non-sense. One should take the cutoff A to be arbitrarily 
large, and make the necessary subtractions of divergent and finite terms. However 
making subtractions is a bit of a pain. Firstly, since the lower order terms in the 
1/c expansion are modified by the inclusion of any higher order term, one must 
recalculate counterterms every time a new higher order term is added. Secondly, 
this means that to restore power counting one must tune coefficients carefully. 

If instead one computes with the multipole expanded lagrangian, the result in 
dimensional regularization is 



iT 



(2) _ c 2 F C{R)g 2 f d D k k 2 



2M 2 d> J (2n)° ( k 2_p +ie ^ E _£_ + ko + ie 



Note that, by design, the whole dependence on c is explicit in the coefficient, 1/c 5 . 
The divergent part is now 

. r(2 ) _ c 2 F C(R)a S rfn ^ (E-p 2 /2M) 3 
div " 2vr [S) M 2 c* 

2.8 NRQCD - An Effective Field Theory: Is All This Really Necesary? 

Although nobody has bothered to proof the validity of NRQCD as an effective 
theory (in the sense of a factorization theorem, as explained above), we believe this 
to be true. Why? 

The answer is that NRQCD has the right ingredients. That is, one could simply 
ask what terms should be contained in a local lagrangian field theory describing 
gluons and non-relativistic quarks. The answer 113 is the lagrangian we have been 
considering, 

C = ~F%F^ + tfm + ^fD 2 )^ + • Bip H (13) 

The coefficient of the F 2 and ifj>D t if terms have been scaled to one by wave function 
renormalization, as usual. The coefficient of the tjr D 2 i[> term defines the mass M. 
Coefficients of higher dimension operators (cf, cd, . . .) remain undetermined. In a 
non-perturbative treatment of the theory, as for example on the lattice, one can 
treat cf,cd,cs as additional free couplings to be determined by comparing with 
experiment (or with full, continuum QCD, in which case we call the procedure 
"matching"). 

It would appear then that all the careful analysis of the previous sections, and in 
particular the multipole expansion, is really unnecessary. The lagrangian of Eq. (|f) 
describes the same physics as our multipole expanded lagrangian of Eq. ([l2|) . The 
only difference is that in the first case one may not systematically study each order 
of the 1/c expansion. This may not be important for many applications, and in 
such cases the simpler formalism should be employed. But, as we will see, this 
is not always the case. An important example is given in the next chapter. The 
calculation of the decay of P-wave quarkonium exhibits infrared divergences in 
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the matching to the effective theory that can be eliminated properly only if one 
includes the correct set of operators. Counting of powers of 1/c singles out the 
correct operators. Moreover, if some of the 1/c effects were left implicit in the 
lowest order states (or elsewhere) one could in principle have some of the infrared 
cancellation take place against part of the state definition: this would make such a 
cancellation practically intractable. 

2.9 Spin Symmetry 
The Lagrangian 

has a U(2) x U (2) symmetry 

Ipa -» #a/#/3 Xa -» TapXff 

with R^R = T^T = 1. Please note that while the field ip annihilates a heavy quark, 
the field \ creates a heavy anti-quark. The abelian factors U(l) x U(l) correspond 
to separate conservation of heavy quark and anti-quark numbers. The SU(2)'s are 
spin symmetries. They are broken by the spin-flip magnetic moment interaction 
that appears in L c -i, and which is of order c~ 5 / 2 in the multipole expansion: 

L c - 1M p= / d 3 x -2E^^(x, t )a- B{6,t)iP(x, t) + ■ ■ ■ 

None of these observations come as a surprise to the reader who has kept in mind 
the molecular physics analogue of quarkonium. 

In non-relativistic quantum mechanics orbital angular momentum L and spin S 
are separately conserved, and therefore so is the total angular momentum J = L+S. 
States are specified by s Lj, where in the standard spectroscopic notation S takes 
values denoting the dimension of the SU(2) representation (S — 1 for the singlet, 
S = 3 for triplet (adjoint), etc), L = S,P,D, . . . and J = 0, 1, 2, ... In the following 
table we give the possible quark-anti-quark states for the lowest few values of L, 
and the name of the state for the case of charm quarks: 



L = 




'So 

3 s, 


Vc 

J/V 


L = 


1 






'Pi 


h c 




3 Po 


XcO 




3 Pi 


Xcl 




3 P2 


Xc2 



In leading order in the 1 / c expansion states are classified this way, and their prop- 
erties are related by spin symmetry. The names of b-quark states are obtained by 
replacing c — > b, except for the 3 Si which is called "T." 

The spin symmetry is most easily applied by use of the Wigner-Eckart theorem. 
This is best explained through examples. 
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Example 1: S-states 



Construct 2x2 matrix H a ^ representing the four 1 Sq, 3 S\ states. Under spin 
symmetry 

H -> RHT ] . 

Then, for any matrix T and any state X free of heavy quarks (and therefore a singlet 
under spin symmetry), the Wigner-Eckart theorem gives 

(X\ X ^\H)=^ x TrTH. 

Here £x is the reduced matrix element; generally it depends on kinematic variables 
and the parameters of the theory. We have used the same letter H to denote the 
specific quarkonium state on the left hand side and the corresponding matrix on 
the right. Explicitly (for 6-quarkonium) 



0=1 



Here and below <t° = diag(l, 1). 
This implies, for example 



(x\x^\v b ) = 

(X|xW|T(e)) = ^ 



where e*is the polarization vector of the spin-1 T. 



Example 2: P-states 

More interesting is the case of the 12 P-wave states. The 12 = 3 x 2 x 2 states can 
be represented by 

H ma0 , m= 1,2,3 a J = 1,2 

and the components are found by insisting they transform correctly under rotations 
(index m) and spin symmetry (a and 0), combined into appropriate J. Without 
proper normalization: 

H maP - O + e mjk al x k i + o^xT + °l^ m 

Here X™ 3 is symmetric and traceless. 

Exercise: (a) Normalize these properly. To this end compute TiH^H and verify 
that the coefficients of all states are unity. 

(b)In the rest frame of the quarkonium state (p = 0) what is (0|x'P^>|if)? Calculate 
(0\x^ D m Tip\H) . Find the relations between these matrix elements of x.J and h. 
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3 Annihilation Decays of Quarkonium 

3.1 Color Singlet Model 

Consider the decay of r\ c into light hadrons. To this effect model the state by a pure 
singlet, i] c = (QQ)i, that is, with no (QQ)$ component. Let us represent the QQ 
in the 1 S'o state as in the quark model 

1 f d^n ffi f mn 

l "^7!^/(^* (,1 ^7^ |c "" (,1^ '" ( ^ 1> • 

where i, j are color indices, m, n are spin indices and ip{q) is the Fourier transform 
of the coordinate space wavefunction 

By spherical symmetry 

1>W = 

where r = \x\. 

We stop here for a small digression. We note that something is missing. Even 
in the molecular physics analogy of Sect. 1 the state would be represented by a 
wave-function <f> n (x,X)i[)(X), with (j> n {x,X) the wavefunction for the fast degrees 
of freedom. We are missing the wavefunction for the glue! We will look the other 
way and carry on, as if the glue component were trivial. 



k,A 



q-kV 



+ crossed + higer order 



WQQ^ B 



Fig. 13. Feynman diagrams for the cc — ► gg amplitude, used in the 
calculation of the rj c width. 



In order to calculate the total inclusive width T(ri c — > hads) we use duality. This 
is the physically plausible assumption that the sum over all final hadronic states 
should be well approximated by the rate into partons, 

T( Vc -» hads) - Yl X ) ~ 99)- 

x 
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With our representation of the r\ c state we then have 
1 f d 3 k 2nS(M ri - 2\k\) 



T(j] c — > hads) 



2M„ J (27r) 3 2fc 



M„ 



Here M. is the invariant (T-matrix) amplitude. Suppressing color and spin indices 

1 f d 3 q 



M( Vc ^g{k)g(-k)) 



:iP(q)M(c(q)c(-q) -> g(k)g(-k)). 



In perturbation theory the amplitude M(c(q)c(— q) — + g(k)g(—k)) is computed from 
the diagram in Fig. O, and is given by 



M(c(q)c(-Q) - 9(k)g(-k)) = -igi(T«T A )v 7 ■ e*(k') 



7 • (q — k) — m, 



■7 • e*(k)u 



where q — (m c + E,q), q' = (m c + E,—q) and k = ^M v (l,h) (\h\ = 1). Note 
that, for now, I have suppressed powers of c which can be restored by dimensional 
analysis. The denominator 



(q - k) 2 -m 2 c Ri -2q ■ k = -M v [(m c + E) - q '■ h] 



(14) 



is dominated by the mass term for \q\ < m c and there it is approximately constant 
(or "flat") and « —2m 2 . The wavefunction ip(q) has support within one inverse 
Bohr radius, 1/aBohr ~ ot s m Cl which is indeed smaller than m c if a perturbative 
value for a s can be used. 



WW 




Fig. 14. The amplitude for cc — > gg is effectively local: when the external 
charm quarks are non-relativistic the internal quark must be off shell by 
an amount of order of the charm mass. 



We therefore approximate the integral by replacing the propagator by its flat 
value, 



T(r] c — > hads) 



1 



(2M„) 



<1 



(2tt) 3 
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d 3 k 2n8(M„ -2|fc|), , y , , , , . , Al , 

" 1 l ±\M(c(0)c(0) ^ g(k)g(-k)\< 



(2n) 3 2k M v 



Exercise: Complete the calculation. Show this gives 



Sit 2 

r(^hads)«^|i?(0)| 2 . 

3.2 i] c Decay in NRQCD 

We would like to carry out the calculation of the previous section in a more rigorous 
manner. To this end we would like to use the effective theory, NRQCD. However, 
several problems immediately come to mind: 

• The effective theory cannot describe properly gluons or quarks with E, \p\ ~ 
m c (because the 1 /c expansion will result in an expansion in velocity v with 
v/c ~ \p\/m c c or v 2 ~ E/m c c 2 ). 

• The effective theory has separate conservation of Q and Q numbers. How can 
it possibly describe QQ annihilation in rj c decay? 




P k+p-p' P 

Fig. 15. Feynman diagram for the forward cc scattering amplitude. 



Fortunately there is a common solution to these problems .EZl'Ea Consider the 
Feynman diagram for cc annihilation in Fig. |l3|. The external quark lines, rep- 
resenting the quark- antiquark pair in charmonium, arc almost on-shell. But the 
internal quark line is off the mass shell by the mass of the quark itself! Effectively 
the interaction is local; see Fig. |l4|. Formally the denominator in Eq. (fL4|) can be 
expanded 

1 1 E-q-h 



m c + E — q ■ h i 

and the leading term gives a local non-derivative quark-gluon interaction for cc 
99- 
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= 21m 



Fig. 16. The optical theorem relates the total cc annihilation cross sec- 
tion, which enters the calculation of the width of the 77c, to the cc for- 
ward scattering amplitude, which can be modeled by a local interaction 
in NRQCD. 



Notice however that the external gluons have energy and momenta of order 
m c . This, of course, is not appropriate for treatment in the effective theory. So we 
are not out of the woods yet. The proposed solution is to use the optical theorem 
which gives the rate for cc — > gg from the imaginary part of the forward cc scattering 
amplitude. The lowest order Feynman diagram for cc — » cc is shown is Fig. |l5|. The 
imaginary part has the internal lines off the mass shell by about m c , so the whole 
loop is local on the scale of \q\ ~ ^Bohr ~ TO cCt s ~ m c v and we may replace it by 
a local interaction. This is shown in Fig. [l7| where k is a constant and the ellipsis 
indicate terms of order \q\/m c . The constant k can be easily computed in the limit 
q = 0. 




Fig. 17. The cc forward scattering amplitude can be modeled by a local 
interaction in NRQCD. The constant k is determined by "matching" to 
QCD. 



To get some better understanding of how this works, let's look at the propaga- 
tors in the box diagram of Fig. ^5l Let k stand for the loop momentum and p and 
p' the external momenta. The quark propagators give 

1 1 



k 2 + 2p ■ k k 2 - 2p ■ k 



2G 



and the gluon propagators 



k 2 (k+p-p') 2 

Now, the imaginary part of the diagram gives a cut on the gluon lines, 

Im , 2 1 ■ = -n5(k 2 ) 
k z + le 

so the product of internal propagators reduces to 

This is easy to interpret: the two gluons are on shell, the two quark propagators 
are just as before, Eq. (|T^). 

In NRQCD the annihilation process is modeled by adding to the lagrangian 
a non-hermitian term that gives the imaginary part of the QQ forward scattering 
amplitude. That is, the lagrangian is augmented by a local four-fermion operator 

with / a dimensionless c constant chosen judiciously to give the right answer to this 
order in the 1/c expansion. This operator is the product of two ( 1 S'o)i bilinears 
(the subindex refers to the color channel, the spin quantum numbers are inside the 
parenthesis). The first bilinear annihilates the ( 1 5o)i cc, while the second creates 
it from the vacuum. 

The constant / is determined by "perturbative matching" . This means that 
the rate should correspond to that of full QCD to the accuracy of the order of 
relevance of the NRQCD approximation. But the same constant / should reproduce 
the perturbative QCD calculation of a partonic cross section. Let's consider this 
procedure in some detail. In QCD the forward scattering amplitude is given by the 
box diagram of Fig. [l5[ The spinors come in as 

^7 M [7 • (p + k) + m]^ v uu^ v [^ ■ (p — k) + m^^v. 

In the NR-limit we can replace 



i and d —— ■ 

0/ y-x 1 

The product of Dirac matrices is reduced by -ff i -f X j L ' = jf^ri + ^^x^ _ ^x^nv 
ie^^P'ypjs and 7^7^ = rf- v — icr^ v . The needed binomials are 



,0 ^ = 

y \'a tp ji = 1 

xH v = 

fi = i 
vcr 0l u = — ixVV 
va n u = 



To count dimensions, recall S = J dtd?xC is dimensionless, C = ifid/dt if), so [ip^] ~ L 
Since % = 1, [Mc] ~ L _1 so we have [ J dt Sxj^-^^ X xH\ ~ L 4 L 2 L~ 6 ~ 1. 
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Notice that when we construct singlet combinations we have, in addition, 

£ x + *ty = o. 

singlet 

Therefore the only non-vanishing singlet-singlet contribution is from 
ie^°ie^ (p + k) x {p - kYxH^X 

The coefficient is just [(p + k) ■ (p — k) — (p + k)°(p — k)°] and this, on-shell (when 
we take the imaginary part of the amplitude), is m\. 

Exercise: Complete the calculation. Show that the spin singlet, color singlet 
part of the cut diagram is 

In NRQCD this is to be reproduced by a term 

j 

— oXV^X- 
mj 

in the Lagrangian. The symbols \ an d "0 are fields, but the tree level amplitude 
corresponds to replacing them by spinors. The matching condition is 

Im f = - 

J 9 

Finally, we compute the rate r(ry c ) in NRQCD. Treating £ as a perturbation, 
the energy shift is d 

AE = / fa# t XX t '0hc) 

r ' c m 2 c c 2M Vc 

The imaginary part of AE is just ^T, so we have finally 

We can find the relation between this result and that of the color-singlet model 
by using the vacuum insertion approximation 

(Vc\^XX^\Vc) « MxH\r)c)\ 2 

with 

(oixVi»?c>»v^ryj«(o) 

We close this section with a short digression. The reader may be concerned with 
the appearance of complex coefficients of hermitian operators in the lagrangian. 

tn the standard formula AE = ^ffij^ g° to relativistic normalization: ('tl>\ip) = 2M^ e <5 3 (0) = 
2M Vc V. The factor of V' 1 is used in H'V' 1 = W = -AC. 
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This renders the Hamiltonian non-hermitian and jeopardizes conservation of prob- 
ability. A little thought shows that this is as expected. The situation is similar to 
that encountered in the effective two level Hamiltonian for K° — K° mixing. There 
the Hamiltonian is given in terms of hcrmitian matrices M and r as H = M + |I\ 
The loss of hermiticity encoded in T arises because the states into which the Ja- 
mesons may decay are not incorporated in the theory. Similarly in NRQCD we 
have truncated the theory by excluding energetic gluons and light quarks, but these 
are the decay products of heavy quarkonia. 



1M» 



Fig. 18. Feynman diagram for cc annihilation entering the calculation of 
the T) c decay width. 
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3.3 Radiative Corrections and Decays of P-Waveonium 

In the cc annihilation diagram of Fig. [l^, used to describe r\ c decay, the giuons 
have momentum |fc| = -M^. Consider now the radiative corrections shown in 



the diagrams of Fig. 19. We are particularly interested in the kinematics were the 
momentum s is soft, that is, small. In both diagrams the sub-diagram corresponding 
to cc — > gg of Fig. |l8| is much as before. That is, the internal quark is far off the 
mass shell, and the two external giuons are hard. As before we attempt to replace 
the interaction by a local ccgg vertex. 




Fig. 19. Feynman diagrams representing soft gluon radiation as a cor- 
rection to the amplitude for cc annihilation of Fig. |l8[ The soft gluon 
carries momentum s. 



Let us analyze the emission of the soft gluon. In the first diagram of Fig. [Jij 
the quark line emitting the gluon gives 

7- (±P + g-s) + m c ^ 1 -{\P + q)+m c 

(±P + q- s) 2 -ml 7 (±P + g) 2 -m 2 

, (1 + 7°) „ (l + 7°)^c 
~ -2s° + s 2 /m c 1 (|P + qy - m 2 c 

1 + 7° \ m c 

2 J (|P + g)2_ m 2_ ' 

We have left the propagator on the right side unexpanded since it will be amputated 
when we compute an amplitude. Similarly, the second diagram in Fig. [l^ gives 

7- (-\P + q) + m c ^ 1 -{-\P + q + s)+m c 
(-\P + q) 2 ~ml 7 {-\P + q + S y-ml 
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( 



1-7' 
2 



o 



) 



+7' 



o 




We see immediately that although each diagram is separately infrared divergent, 
ie, they contain a singularity as s° — > 0, the divergence cancels in the sum. Phys- 
ically the cancellation can be easily understood: the coupling involves 7 , which 
corresponds to the "charge" of the quark, but c and c have opposite charges. 

To see that this physical intuition is correct, let's compute the next term, in- 
volving the current 7*. Using |P + q = (m c c, 0) + {\E, q) and neglecting q° and s, 
the Dirac structure of the soft gluon-quark interaction is 



The leading term (for fj, = 0) is what we computed above. The new terms of interest 
are the ones of order m c . The second diagram gives, similarly 



So we see that there is at this order an infrared divergence in these diagrams. It 
appears only in the y, = i terms. However, it does not afflict decays of S-wave 
charmonium since the factor q l changes L by one. e Formally, the operator with the 

divergent coefficient is Dip, rather than x V*- But 



[7 ' ( 2 P + 1 - s ) + m c] 7 M [7 ' ( 2 P + l) + m c] 



[m c (l + 7 °) +7 • q] 7 M [m c (l + 7 ) + 7 • q 




[7 ' (- 2 P + 9) + mc] Y [7 • (- + <? + s ) + m c\ 




(Olx^l^c) = 0. 



However L = 1 operators do interpolate for P-wave quarkonium. 



Exercise: Using the methods of Sec. 2.£ show that 



(0\ X ^-D^\h c ) =0 
(0|x^-^|xi} = 



and find the relation between (0|x^cr • Dip\xo) and (0|x^<? ■ D^|x2)- 



The amplitude has q ■ e, where e is the gluon polarization vector. 
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In the color singlet model there seems to be no way out of this problem with 
infrared divergences in the decay of P-wave quarkonium. Let us emphasize that 
the problem is that we have found an infrared divergence in the calculation of the 
matching of a coefficient, which should be given by short distance physics. Usually 
one finds a way to shove the infrared divergences into matrix elements of operators, 
but there is no candidate for an operator into which to shove this infrareiLdivergence. 

However, in NRQCD there is a neat resolution to this problemEa In leading 
order, the decay is described by local interactions 

sc = Ii^o 1 ( 3 p ) + M^ilo 8( 3 5l) (15) 

where f's are coefficients determined by matching and O are operators. The sub- 
script denotes the relative color of the quark bilinears, singlet or octet, and in 
parenthesis we have their J L$ number. The operators are 

8 ( 3 S 1 )=^(a i T a )xxH<r i T a )iP 




Fig. 20. Feynman diagram for cc annihilation into light quarks. It rep- 
resents the leading contribution to the matching of the color octet oper- 
ator C 8 ( 3 Si). 



Exercise: Show that at tree level the diagram of Fig. |18| gives the matching to 
the P wave operator 

Im/l(3po) ^^T a ' (16) 

while the diagram of Fig. [2^, where there are rif species of light quarks in the 
possible final state on the right, gives the iS-wave coefficient 

Im/ 8 ( 3 ^) = ^ (17) 

3-4 Power Counting 

Decays of 5-wave Quarkonium 

For this discussion to make any sense we must show that the two operators in 
Eq. JIe| ) contribute to the same (leading) order in 1/c. It is instructive though to 
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begin by revisiting the decay of the S'-wave quarkonium. Let us set the normal- 
ization arbitrarily by saying that the contribution from the operator Oi( 1 5'o) is of 
order one: 

OiCSo): (Vc\tfxxH\ric)~l 
Consider other possibilities for 8C: 

Oi( 3 Si) : (*7<#V X xVV|r?c) 

This vanishes by spin-symmetry. However, we can include, perturbatively, interac- 
tions that break spin-symmetry and then the result need not vanish. Recall 

Now a ■ B flips spins, but also changes relative color from singlet to octet. One 
can restore the quarks into a singlet color combination by inserting, in addition, an 
electric dipole (El) operator 

Thus we estimate 
Consider next 

OsCS ) : (Vc\1> 1 T a xx ] T a 1>\Vc)- 

Inserting the El operator twice (once is not enough because it has L = 1), we 
estimate its order to be 

'.9 2 V (9 2 \ 1 



c 3 J \ c J c 4 

The operator C 8 ( 1 Pi) requires only one El insertion, but comes built in with a 1/c 2 
suppression: 

For these non-perturbative matrix elements ahe strong coupling constant should 
be taken large, g 2 /c <~ a s ~ 1. We conclude that the corrections to the T] c decay 
width from C 8 ( 1 S'o) and O^Pi) are suppressed by 1/c 4 , and from C 1 ( 3 S' 1 ) by 
1/c 6 . The leading correction is in fact of order 1/c 2 and arises from an operator 
with the same quantum numbers as the leading operator, but with two derivatives: 



("he matrix elements arc considered to all orders in the zeroth order lagrangian of NRQCD. 
We have used perturbation theory only with respect to the dimensionful parameter 1/c. 
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Decays of P-wave Quarkonium 

We now perform an analogous analysis for the case of decays of P-wave states. The 
operators we focus on are the ones in the effective lagrangian of Eq. (p"5|). Oi( 3 P ) 
has just the right quantum numbers to interpolate the color singlet P-wave states, 

-L-^Io^p^p) 1 



c 3 

while Os( 3 Si) requires an El transition, 

m z c c Vr' / \ c / e 3 

We see that these two operators contribute to the P-wave decay width at the same 
order in the 1/c expansion. It is easy to see that there are no other operators that 
contribute at this order. 

We can now use the Wigner-Eckart theorem to express the decay widths in 
terms of reduced matrix elements 

ln \ KQIx^-^lxco)! 2 

{<Dl) = 2M~ X 

{ ° 8) = 2M X 



Using the results of matching at leading order in Eqs. (|16|) and ( |l7|) the widths for 
the four P wave particles are given in terms of these two reduced matrix elements: 



6m 



TO 4 



3m? 



6m A c 

We are now in a position to understand the cancellation of infrared divergences. 
In the expression 

Imh( s P ) //n x , Im/ 8 ( 3 5i) /<n , 
r Uco) = 7-3 — Ci H 2 

the infrared divergence in Im/i( 3 Po) is precisely cancelled by a divergence in the 
matrix element (Os)- The resolution of this problem is responsible for much of the 
renewed interest in NRQCD. 
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3.5 Decays: Summary 



In computing decays of quarkonium into light hadrons we include in the effective 
lagrangian local four-fermion operators 

S£ = Y,fnO n 

n 

and compute 

ir = ^im/„(o„), 

n 

where Im/ n is computed by matching to QCD. The relative importance of these 
terms is determined by power counting in 1 /c, using the minimal number of required 
interactions from the multipole expanded NRQCD effective lagrangian, and using 
g 2 /c~ a - 1. 

The operators O n are of the form 

«— > 

with n m a tensor product of D, T a and a 1 . 



4 Concluding Remarks 

There are many open-questions that remain un-addressed. In the standard presen- 
tation of the subject 113 the color octet bilinears in the lagrangian of NRQCD, such 
as Os( 3 Si) of Eq. (|l5|), are interpreted as literally producing and annihilating color 
octet states. It is not quite clear what these states are. It is reasonable to assume 



they correspond to the hybrid states p£ Sec. 1.4. But doing so would invalidate the 
standard analysis for power countingjla which requires the size of this states to still 
be of the order of asohr rather than the bag radius Aq^ d . 

A closely related question is whether the non-perturbative potential, which 
grows linearly with distance at large distances, should be used in NRQCD. The 
Born-Oppenheimer approximation certainly suggests this is the right approach. 
However, it remains possible that the condensates of the higher dimension opera- 
tors fully incorporate the non-perturbative information relevant to all the questions 
that NRQCD attempts to answer. In that case, including both the non-perturbative 
potential and the condensates of higher dimension operators would result in double- 
counting. 

Another question that remains unsolved is the rather arbitrary assumption 
that g 2 /c must be taken as of order unity in order for power counting to work, as 



presented in Sec. 3.4. It is not enough to insist that g 2 ~ c in general, for this must 
not be the case in calculations of matching coefficients of the operators in 8C. The 
special rule g 2 /c ~ 1 applies only when calculating matrix elements of operators. 

Because of time and space limitations I have not included in these lectures a 
discussion of the analysis of quarkonium production in NRQCD. Yet most of the 
ongoing work on the, subject is in this realm. There are many good reviews the 
reader may consult S3 
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